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ABSTRACT 

High-frequency twin peak quasiperiodic oscillations (QPOs) are observed in four microquasars, i.e., Ga- 
lactic black hole binary systems, with frequency ratio very close to 3 : 2. In the microquasar GRS 1915+105, 
the structure of QPOs exhibits additional frequencies, and more than two frequencies are observed in the 
Galaxy nuclei Sgr A*, or in some extragalactic sources (NGC 4051, MCG-6-30-15 and NGC 5408 X-1). The 
observed QPOs can be explained by a variety of the orbital resonance model versions assuming resonance of 
oscillations with the Keplerian frequency Vk or the vertical epicyclic frequency Ve, and the radial epicyclic 
frequency v,, or some combinations of these frequencies. Generally, different resonances could arise at 
different radii of an accretion disc. However, we have shown that for special values of dimensionless black 
hole spin a strong resonant phenomena could occur when different resonances can be excited at the same 
radius, as cooperative phenomena between the resonances may work in such situations. The special values 
of a are determined for triple frequency ratio sets Vk : Ve : Vr = i' : ( : « with s,t,u being small integers. The 
most promising example of such a special situation arises for black holes with extraordinary resonant spin 
a = 0.983 at the radius r = 2.395 M, where Vk : Ve : v^ = 3 : 2 : 1. We also predict that when combinations 
of the orbital frequencies are allowed, QPOs with four frequency ratio set 4 : 3 : 2 : 1 could be observed in 
the field of black holes with a = 0.866, 0.882 and 0.962. Assuming the extraordinary resonant spin a = 
0.983 in Sgr A*, its QPOs with observed frequency ratio « 3 : 2 : 1 imply the black hole mass in the interval 
4.3 X 10^ Mr: < M < 5.4 X 10^ Mq, in agreement with estimates given by other, independent, observations. 
Key words: Accretion, accretion disks - Black hole physics — X-rays: general 

1 Introduction 

Quasiperiodic oscillations (QPOs) of X-ray brightness had been observed in many 
Galactic low-mass X-ray binaries containing neutron stars (see, e.g., van der Klis 2000, 
2006, Barret et al. 2005, Belloni et al. 2005, 2007) or black holes (see, e.g., McCHntock 
and Remillard 2004, Remillard 2005, Remillard and McClintock 2006). Some of the 
QPOs are in the kHz range and often come in pairs (Vupp,Vdown) of twin peaks (often 
called double peaks) in the Fourier power spectra. Since the peaks of high frequencies 
are close to the orbital frequency of the marginally stable circular orbit representing the 
inner edge of Keplerian discs orbiting black holes (or neutron stars), the strong gravi- 
ty effects have to be relevant in explaining high-frequency QPOs (Abramowicz et al. 
2004). 

The twin peak QPOs were observed in four microquasars, namely GRO 1655-40, 
XTE 1550-564, H 1743-322, GRS 1915H-105 (Torok et al. 2005). In all of the four 
cases, the frequency ratio of the twin peaks is very close to 3 : 2 and the orbital reso- 
nance model assuming non-linear resonances between oscillations in Keplerian Vk and 
epicyclic frequencies Ve (vertical), or Vr (radial), or their combinations, seems to be the 
potential explanation of the microquasar kHz QPOs. The resonances of oscillations 
with the epicyclic frequencies were firstly mentioned in Aliev and Galtsov (1981). In 
the context of QPOs observed in the neutron star and black hole systems, the orbital 
resonance model was introduced and developed by Kluzniak and Abramowicz (2000), 
see Abramowicz et al. (2004). 

According to the resonance hypothesis, the two modes in resonance should have 
eigenfrequencies Vr (equal to the radial epicyclic frequency) and Vv (equal to the ver- 
tical epicyclic frequency Ve or to the Keplerian frequency Vk). While models based 



2 



on the parametric resonance identify the two observed frequencies of the twin peaks 
(Vupp, Vdown) directly with the eigenfrequencies of a resonance, models based on the 
forced resonance allow to observe combinational (beat) frequencies of the modes, or 
oscillations with the combinational frequencies entering the resonance (Landau and 
Lifshitz 1976). Both parametric and forced resonance models make clear and precise 
predictions about the values of observed frequencies in connection with spin and mass 
of the observed object, at least in the case of black holes (Torok et al. 2005). 

In some sources, more than two high-frequency peaks are observed. The micro- 
quasar GRS 1915+105 reveals high-frequency QPOs appearing at four frequencies 
with the lower and upper pairs in the ratio close to 3 : 2 (Remillard and McClintock 
2006) and even a fifth frequency was reported (Belloni et al. 2001). An additional 
sixth frequency was mentioned (Strohmayer 2001), although not confirmed. In Sgr A*, 
three frequencies were reported with ratio close to 3 : 2 : 1 (Aschenbach 2004, Aschen- 
bach et al. 2004, Torok 2005). In the galactic nuclei MCG-6-30-15 and NGC 4051, 
two pairs of QPOs were reported with the ratios close to 3 : 2 and 2:1, respectively (La- 
chowicz et al. 2006). In the source NGC 5408 X-1 oscillations with three frequencies 
of ratio close to 6 : 4 : 3 were observed (Strohmayer et al. 2007). 

In the microquasar GRS1915+105, an almost extreme Kerr black hole with awl 
is expected (McClintock et al. 2006), and all the five (six) frequencies of QPOs can be 
explained in the framework of the extended resonance model with the hump-induced 
oscillations, predicting the black hole spin a = 0.9998 and its mass M = 14.8 M© 
(Stuchlfk et al. 2006, 2007d). In the extended resonance model, forced resonances 
of the epicyclic oscillations with an additional oscillation induced by the humpy or- 
bital velocity profile (related to the physically privileged locally non-rotating frames) 
that occurs in Keplerian discs orbiting Kerr black holes with a > 0.9953 are assumed 
(Aschenbach 2004, 2006, Stuchlfk et al. 2004, 2005, 2007c). In the "humpy" extended 
resonance model, all the oscillations in resonance can occur at, and be related to the 
exclusively defined, "humpy radius" with extremal orbital velocity gradient within the 
humpy profile (Stuchlfk et al. 2007c). However, this model can be relevant only for 
near-extreme Kerr black holes with spin a > 0.9953 for Keplerian discs, and with even 
higher spin a > 0.9998 for marginally stable thick accretion discs (Stuchlfk et al. 2004, 
2005). 

In order to explain complex frequency structures observed in some black hole sys- 
tems, we expect that more than one resonance occur in a Keplerian disc, and different 
versions of the orbital resonance model could be acting. Of course, one version of the 
orbital resonance model appearing at different resonant points, at different radii of the 
disc, is also possible (Stuchlfk and Torok 2005). Note that the so called total preces- 
sion resonance model, assuming resonance of oscillations with the Keplerian frequency 
Vk and the total precession frequency Vt = Ve — Vr, can explain quite well the multi- 
resonant phenomena observed in the neutron star atoll source 4U 1636-53 (Stuchlfk et 
al. 2007e) and in some other atoll sources (Bakala et al. 2008, in preparation). 

Generally, two resonances could occur at different radii, and for special values of 
the dimensionless spin a, the bottom, top, or mixed frequencies coincide, reducing two 
frequency pairs into a triple frequency set (Stuchlfk et al. 2007b). Physically more in- 
teresting seems to be the case, when a triple frequency set Vk, Ve, Vr with rational ratio 
arises at a given radius. The resonant phenomena can be then expected stronger than in 
the case of internally independent resonances occurring at two different radii. Clearly, 
while sharing a given radius, the resonances could be causally related and could co- 
operate efficiently (Landau and Lifshitz 1976). That is the reason why we focus here 
attention on the possibility of causally related resonances appearing at a fixed radius. 
Such a situation is possible only in the field of Kerr black holes with special values of 
the dimensionless spin a, when the value of a is related to the corresponding triple fre- 
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quency ratio set, and concrete versions of the resonance model that are realized. Since 
the strength of the resonance and the resonant frequency width decrease rapidly with 
the order of the resonance n + m (see Landau and Lifshitz 1976), in the following we 
restrict ourselves to the cases with n,m<5. 

2 Orbital Resonance Model 

The standard orbital resonance model (Abramowicz et al. 2004, Torok et al. 2005) as- 
sumes oscillations of an accretion disc orbiting a rotating black hole described by the 
Kerr geometry, or a neutron star that could be represented by the Schwarzschild ge- 
ometry, or, more precisely, by the Hartle-Thorne geometry (Hartle and Thorne 1968). 
The accretion disc can be approximated by a thin disc with Keplerian angular velocity 
profile, or by a thick toroidal disc with angular velocity profile given by distribution 
of the specific angular momentum in the fluid of the toroidal disc. The frequency of 
the disc oscillations is related to the Keplerian frequency (orbital frequency of tori), 
or the radial and vertical epicyclic frequencies of the circular test particle (geodetical) 
motion. The epicyclic frequencies can be relevant both for the thin, Keplerian discs 
with quasicircular geodetical motion (Kato et al. 1998, Novikov and Thorne 1973) and 
for thick, toroidal discs with non-geodetical quasicircular rotation kept by the pressure 
gradients of the tori (Schnittman and Rezzolla 2006, Rezzolla et al. 2003, Rezzolla 
2004ab). However, with thickness of an oscillating toroid growing, the eigenfrequen- 
cies of its radial and vertical oscillations deviate from the epicyclic test particle fre- 
quencies (Sramkova 2005). Here we focus our attention to the Keplerian thin discs. 

Different versions of the orbital resonance model could be classified according to 
the following criteria: 

a) the type of the resonance (parametric or forced), 

b) the type of oscillations entering the resonance, 

c) the presence of beat, combinational frequencies. 

Thus, according to the first criterion, two main groups of orbital resonance model ver- 
sions exist, differing by the type of the resonance. In both of them, the epicycUc fre- 
quencies of the equatorial test particle circular motion have a crucial role (Torok et al. 
2005). 

The internal resonance model is based on the idea of parametric resonance between 
vertical and radial epicyclic oscillations with the frequencies Ve = coe /27t and Vr = 
(Or/27t. The parametric resonance is described by the Mathieu equation (Landau and 
Lifshitz 1976) 

5e-i-(»^ [i-F/jcos((OrO]8e=o. (1) 

Theory behind the Mathieu equation implies that a parametric resonance is excited 
when 

COr Vr 2 

7^ = - = -, «=1,2,3,... (2) 

and is strongest for the smallest possible value of n (Landau and Lifshitz 1976). Be- 
cause Vr < Ve near black holes, the smallest possible value for the parametric resonance 
is « = 3, which means that 2ve = 3Vr. This explains why the 3 : 2 ratio is commonly 
observed in the black hole systems, assuming Vupp = Ve and Vdown =Vr. Note that for the 
internal resonance the oscillating system conserves energy (Landau and Lifshitz 1976). 
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Versions of the resonance model based on the forced resonance come from the idea 
of a forced non-linear oscillator, when the relation of the latitudinal (vertical) and radial 
oscillations is given by the formulae 

56 + (Oq56 + [non-linear terms in 56] =g{r) cos((Oo t) , (3) 

5r + (0^dr+ [non-linear terms in 50, 5r] =/z(r)cos((Oo/)) (4) 

with 

©e = I - ) cOr, (5) 



..1, 

where p,q are small natural numbers and ©o is the frequency of the external force, e.g., 
the gravitational perturbative forces are discussed, for the case of a neutron star with 
"mountains" or accretion columns, and a binary partner of the neutron star or a black 
hole, in Stuchhlc et al. (2008, 2007a), Stuchlik and Hledik (2005). The non-linear 
terms allow the presence of combination (beat) frequencies in resonant solutions for 
50(f) and 5r(f) (see, e.^.. Landau and Lifshitz 1976), which in the simplest case give 

CO- = COe - COr, (0+ = COe + COr. (6) 

Another, so called "Keplerian" resonance model, takes into account possible para- 
metric or forced resonances between oscillations with the radial epicyclic frequency Vr 
or the vertical epicyclic frequency V9, and the Keplerian orbital frequency Vk- 

Such resonances can produce the observable frequencies in the 3 : 2 ratio as well as 
in other rational ratii (note that one of the cases which gives 3 : 2 observed ratio is also 
the "direct" case of p:q = 3:2 corresponding to the same frequencies and radius as 
in the case of 3 : 2 parametric resonance). Therefore, we shall consider both the direct 
and simple combinational resonances. 

The resonance conditions of the parametric and direct forced resonance are com- 
mon, however, physical details, as the time evolution of the resonance, the strength of 
the resonance and the width of the resonant frequencies, are different (see Landau and 
Lifshitz 1976). 

The width of the resonant frequencies in forced resonances differs in the cases of 
direct, sub(super)harmonic and combinational resonances, and depends on the external 
force strength and the damping and non-linear terms. Considering a simple case of ex- 
ternal force of harmonic character with frequency £1 influencing a non-linear oscillator 
with a single degree of freedom and variable u, having eigenfrequency coq, damping 
parameter /i and cubic non-Unear term characterized by parameter a, the equation of 
motion reads (Nayfeh and Mook 1979) 

u + (S)lu = -2ejju-eau^ + ekcos{Q.t) (7) 

where e is a small parameter For primary resonance with Q. ^ coo in the linear regime, 
where both damping and non-Unearities are negligible, the amplitude growing is un- 
bounded according to A ~ ? , but i2 = (Oq exactly (Landau and Lifshitz 1976). The linear 
growing of amplitude is limited by damping and non-linearities and detuning of the 
frequencies is allowed. We describe the frequency detuning by a detuning parameter a 
due to 

i2 = (»o-l-ea. (8) 

The frequency-response equation, relating the amplitude of the resulting oscillations A 
and o in dependence on the amplitude of the excitation k, and the parameters /u and a, 
can be put into the form (Nayfeh and Mook 1979) 

1/2 

■ (9) 



3a ^2^f 2V 
8^^"^ ^[-^ ^ ) 
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The linear response (a = 0) is symmetric, while a > introduces an asymmetry - in 
both cases the peak amplitude 

An,ax=:r^; (10) 

(see Nayfeh and Mook 1979 for details). The maximal frequency detuning allowed for 
resonant phenomena represents few percent of the eigenfrequency coq. 
In the case of superharmonic resonances 

3i2 = (jOo+ea, (11) 

the frequency-response equation reads 



. aA^ 3a . , , / a^A^ 



<J = 3— + ^A^±^-^-/ , (12) 



coo 8coo V»o'^^ 

where 

k 



1/2 



2((»2-n2)' 

The peak amplitude of the resonant overtone oscillation is 

_ aA^ 

COo/i 



(13) 



(14) 



i.e., it depends on the magnitude of the non-linear term, contrary to the case of primary 
resonances. Similar frequency-response equations can be given for the subharmonic 
resonance (il = 3coo + eo) or for combinational resonances (Nayfeh and Mook 1979). 
Generally, the maximal frequency scatter represents few percent of COq again (Nayfeh 
and Mook 1979). Similarly, in the case of the quadratic non-linearity (— eaw^) intro- 
ducing the sub(super)harmonic frequencies Q. ~ 2a)o (ii ^ 5W0). 

In the case of the parametric resonance, where the excitation appears as a time- 
dependent coefficient in oscillatory equations, and for internal resonances between cou- 
pled oscillatory modes, the frequency detuning is given in a more complex way (Nayfeh 
and Mook 1979, Landau and Lifshitz 1976) but the maximal relative frequency detun- 
ing is again in percents and sharply decreases with the order of the resonance (Landau 
and Lifshitz 1976). We shall not go into details leaving them to future work. Here we 
concentrate on the resonant conditions. 

The formulae for the vertical epicycUc frequency Ve and the radial epicycUc fre- 
quency Vr take in the gravitational field of a rotating Kerr black hole (with the mass M 
and dimensionless spina) the form (e.g., Aliev and Galtsov 1981, Kato a/. 1998) 

v^ = aev|, v? = a,v| (15) 

where the Keplerian frequency and related dimensionless epicyclic frequencies are 
given by the formulae 



Vk = 



1 /GM^ '/^ 

27t 



ae = l-4ax ^^^ + 3a'^x ^, 

a, = l-6x-^ + Sax-^/^-3a^x-^. (16) 



Herex=r/(GM/c^) is the dimensionless radius, expressed in terms of the gravitational 
radius of the black hole. 
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For a particular resonance n : m the equation 

nvr = mvv; Vv e {ve, Vk} (17) 

determines the dimensionless resonance radius x„:,„ as a function of the dimensionless 
spin a in the case of direct resonances that can be easily extended to the resonances with 
combinational frequencies. From the known mass of the central black hole (e.§., low- 
mass in the case of binary systems or high-mass in the case of supermassive black 
holes), and the observed twin peak frequencies (Vupp, Vdown), the Eqs (fTSl) - (fTTIi im- 
ply the black hole spin a for various versions of the resonance model, with the beat 
frequencies taken into account (Torok et al. 2005). 



3 Black Holes Allowing Multiple Resonances 

Generally, the resonances could be excited at different radii of the accretion disc under 
different internal conditions; such a situation is discussed in detail by Stuchlfk et al. 
(2007b). On the other hand, physically more interesting case arises when two (or 
more) resonances are excited at a common radius because of probable internal physical 
connection of those resonances. Of course, a causal connection of two resonances 
could hardly be relevant, but is not excluded, if they appear at distant radii. 

A simple possibility of two resonances sharing a common radius can be conceived 
supposing a direct resonance and some combinational resonances of the frequencies 
entering the direct resonances. For example, in the field of any Kerr black hole a direct 
resonance of the epicyclic frequencies Ve : Vr = 3 : 2 admits combinational resonances 
Ve : (ve - Vi) = 3 : 1 and Vr : (ve - Vi) = 2 : 1 (see, e.g., Stuchlfk et al. 2007b). However, 
this kind of behaviour is not exhibited by any of the four microquasars investigated at 
present time (Torok et al. 2005). 

Nevertheless, it is important that there exists a possibihty of direct resonances of 
oscillations with all of the three orbital frequencies, characterized by a triple frequency 
ratio set 

VK:Ve:Vr = .?:f :m (18) 

with s > t > u being small integers. The frequency set ratio (Eq. [TS] ) can be realized 
only for special values of the black hole spin a. The black hole mass is then related to 
the magnitude of the frequencies. 

In order to look for the special values of the black hole spin a and the related fre- 
quency ratios (Eq.fTsll. we have to summarize properties of the Keplerian and epicyclic 
frequencies of the circular geodesic motion in the Kerr spacetimes. It is well known 
that Vk and Ve are defined down to the photon circular geodesic located at Xph that is 
determined by the condition 

fl = flph(x) = ^(3-x) (19) 

while Vr is defined down to the innermost (co-rotating) stable circular geodesic at jCms 
given by 

a = a^,{x)=^(4-V3^y (20) 

Clearly, we have to search for the integer ratios (Eq. [TS] ) at x >Xms. 

The Keplerian frequency VK('",fl) is a monotonically decreasing function of the ra- 
dial coordinate for any value of the black hole spin. On the other hand, the radial 
epicyclic frequency has the global maximum for any Kerr black hole. The vertical 
epicyclic frequency is not monotonic, if the spin is sufficiently high (see, e.g.,Kato 
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et al. 1998, Perez et al. 1997). For the Kerr black hole spacetimes, the locations 
r(a), e(fl) of maxima of the epicyclic frequencies Vr, Ve are implicitly given by 
the conditions (Torok and Stuchlik 2005a) 



P,(x,fl) 

Pr(x,fl) 



1 



2x^1^ + a 
1 2a 

1. 



CLj [x, a) where j G {r, 0} , 

2 



X- 
a 



a 

2 



a 



(21) 



For any black hole spin, the extrema of the radial epicyclic frequency %^ r(fl) must 
be located above the marginally stable orbit. On the other hand, the latitudinal extrema 
e(«) are located above the photon (marginally bound or marginally stable) circular 
orbit only if the limits on the black hole spin a > 0.748 (0.852, 0.952) are satisfied 
(Torok and Stuchlik 2005b). In the Keplerian discs, with the inner boundary x\a^x^^, 
the limiting value a — 0.952 is relevant. 

From the point of view of the observational consequences, it is important to know, 
for which frequency ratios n : m the resonant frequencies Ve(a,n:m), considered as a 
function of the black hole spin a for a given frequency ratio n : m, has a non-monotonic 
character. A detailed analysis (Torok and Stuchlik 2005a) shows that Ve(fl,n:m) has 
a local maximum for n:m> 11:5, i.e., in physically relevant situations (n, m small 
enough for the resonance), it occurs for the ratios Ve : Vr = 3: 1, 4: 1, 5:2, 5: 1. 

Assuming two resonances Vk : Ve = : f and Vk : Vr = i : m occurring at the same x, 
we arrive to the conditions 



ae(fl,x) 



ar(fl,x) 



(22) 



(23) 



that have to be solved simultaneously for x and a. The solution is given by the condition 



a^{x,t /s) —a^{x,u/s) 



where 



a^{x,t /s) 



(24) 



(25) 



(26) 



It is possible to find an explicit solution determining the relevant radius for any triple 
frequency set ratio s:t:u 




x(i,f,M) 



6^^ 



6*2 ± 2V2 V(r -«)(? + u) (3^2 - f 2 _ 2m2) -{f- + 5u^)' 



(27) 



Clearly, the condition f ^ -|- 2u^ < 3s^ is always satisfied. The corresponding black 
hole spin a is then determined, e.g., by Eq. ( |25l l giving a^{x{s,t,u),t /s). Of course, 
we consider only the black hole cases when a < 1 . This condition puts a restriction on 
allowed values of i,f , m. 

The solutions have been found for frequency ratios with s <5. We were looking 
for the black hole spin a allowing ratios s:f:M = 3:2:l, 4:3:2, 4:3:1, 4:2:1, 5:4:3, 
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Fig. 1 . Black holes enabling interrelated strong resonances. The special cases of the triple frequency ratio sets 
of the orbital frequencies Vk, Vg, v, (black solid lines) with the corresponding spin a and the shared resonance 
radius. For coinpleteness we present the relevant simple combinational frequencies Ve — Vr, Ve +Vr, Vk — V9, 
Vk — Vr (grey dashed lines). Notice that the extraordinary resonant spin 03:2:1 =0.983 represents the only 
case when the combinational and direct orbital frequencies coincide at the shared resonance radius. 



5:4:2, 5:4:1, 5:3:2, 5:3:1, 5:2:1. Higher values of the s,t,u are presented and 
discussed in a more general context of triple frequency set ratios that could arise when 
two resonances occur at two different radii, but for special values of the black hole spin 
two of the frequency levels (top, bottom, mixed) coincide (see Stuchlfk et al. 2007b). 

We have shown that the direct resonances could result only in the triple frequency 
sets 3:2:1, 4:3:1, 5:4:2, 5:4: 1, 5:3:1 with values of the spin and the radius Xyj-^ given 
in Fig. [T^-e. The other ratios can be realized with the combinational frequencies 
involved in the triple frequency sets 4:3:2, 4:2:1, 5:4:3, 5:3:2, 5:2:1 as shown in 
Fig. [if for the case of the set 4:3:2. We can summarize the results in the following way. 

a) VK:V9:Vr = 3:2: 1 (see Fig.[T]i) 

fl3;2;l =0.983, Jt:3:2:l = 2.395, Jt:ms = 1.571. (28) 

This case of the so called extraordinary resonant spin ay,2:i =0.983 may allow 
strong resonances, in situations when the Keplerian and epicyclic frequencies are 
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in the lowest possible ratio at the common radius X3:2:i =2.395 (see Fig.[T^), and 
is thus of special interest. In fact, this case involves the lowest structure of direct 
resonances with Vk : Vr = 3 : 1, Vk : Ve = 3 : 2, Ve : Vi- = 2 : 1. Notice that in this 
special case also any of the simple combinational frequencies coincides with one 
of the frequencies Vk, Ve, Vr, and are in the fixed small integer ratios 

Vk 3 Vk 3 Ve 2 ^^^^ 



Ve-Vr r Vk-Vi- 2' Ve-Vi- 1' 

It should be stressed that this is the only case when the combinational frequencies 
not exceeding the Keplerian frequency are in the same ratios as the orbital fre- 
quencies. We obtain the strongest possible resonances when the beat frequencies 
enter the resonances satisfying the conditions 

^ = ^ = 1, ^^ = ^- = U ^^ = 1. ^=.1. (30) 
Vk 3 Vk-Vi- 2 VK-ve Vr 

b) VK:ve:Vr = 4:3: 1 (seeFig.IIh) 

04.3.1=0.866, JC4:3:1=2.88, =2.539. (31) 

In this case the combinational frequencies give additional frequency ratios 

Vk 4 Ve 3 Vk 4 2 Ve-Vr 2 



Vk - Vr 3 ' Ve - Vr 2 ' Ve - Vr 2 1 ' Vr r 
and the strongest resonant ratios 



(32) 



^ = ^ = 1, ^ = ^ = 1, ^^ = 1. (33) 

Vk 4 VK-Vr 3 VK-Ve 

Here, using the combinational frequency Ve — Vr, we obtain the other three fre- 
quency ratio sets (see Fig. [TJ) 

VK:ve:(ve-Vr)=4:3:2 (34) 

and 

VK:(ve-Vr):Vr=4:2:l. (35) 

Now, the four observable frequency ratio set is possible when the combinational 
frequency is mixed with the orbital frequencies 

(vk = Ve + Vr) : (ve = Vk - Vr) : (ve - Vr) : (Vr = Vk - Ve) = 4 : 3 : 2 : 1 . (36) 

c) Vk: Ve :Vr = 5 : 3 : 1 (see Fig.[ll;) 

a5.3.i =0.962, X5:3:i =2.083, xms = 1.820, (37) 
and the combinational frequencies (not exceeding Vk) are in the ratios 
Vk Vk 5 VK-Vr 4 Ve 3 



VK-Vr Ve+Vr 4' Ve 3' Ve-Vr 2' 

VK-Ve _ Ve-Vr _2 Vk 5 

Vr Vr 1 ' Ve-Vr 2 ' 



(38) 



Then we can generate triple frequency sets involving the combinational frequen- 
cies with the ratios in the form 



VK:(ve-Vr):Vr = 5:2:1 



(39) 
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and 

VK:V9:(ve-Vr)=5:3:2. (40) 

Moreover, using the combinational frequencies we could obtain two sets of four 
frequency ratios 

Vk: (VK-Vr=ve+Vr):ve:vr = 5:4:3: 1, 

VK:Ve: (ve-Vr=VK-Ve):Vr = 5;3;2; 1, (41) 

and one set of five frequency ratio 

Vk: (VK-Vr = ve+Vr):ve: (v9-Vr = VK-Ve):Vr = 5:4:3:2: 1. (42) 

d) VK:Ve:Vr = 5:4:l (see Fig.[I]l) 

a5;4;i =0.775, JC5:4:i =3.240, .j:ms = 3.033. (43) 
The combinational frequencies can give the ratios 

Vk 5 Ve 4 Vk 5 Ve-Vr 3 



Vk - Vr 4 ' Ve - Vr 3 ' Ve - Vr 3 ' Vr 1 

and the strongest resonant ratios 



Vk - Vr 3 ' Vk - Vr 3 ' Ve-Vr 2 r 

Vr _ 2 Vk 5 

VK-Ve r Ve-Vr 2' 

and the strongest resonant ratio 



Vr 2 

This case leads to the triple frequency ratio sets 



(44) 



Xi±^ = ^ = l, ^^^ = 1 = 1, ^^^- = 1. (45) 

Vk 5 VK-Vr 4 VK-Ve 

This case leads to the triple frequency ratio sets 

VK:(ve-Vr):Vr = 5:3:1, 

VK:Ve:(ve-Vr)=5:4:3. (46) 
Here, only one of the four frequency ratio sets is possible, namely 

(vk = Ve + Vr) : (ve = Vk - Vr) : (ve - Vr) : Vr = 5 : 4 : 3 : 1 . (47) 

e) VK:ve:Vr = 5:4:2(seeFig.[Tt) 

a5.4.2= 0.882, X5:4:2 = 3.407, =2.438. (48) 

Here, the combinational frequencies give the ratios 

Ve 4 Vk 5 Ve 4 2 



(49) 



= ? = 1. (50) 



Ve:(vK-Vr):Vr=4:3:2, 

VK:ve:(VK-Vr)=5:4:3, (51) 
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Fig. 2. Calculation of the width of the resonant radius where the strong resonant phenomena are possible. 
Black lines represent the ratio n:m = VK:Vr, Ve:Vr, dashed lines illustrate the relative frequency detuning 
5 = 5% (i.e., fi:m = 2±0.05 and 3±0.05). Gray solid lines represent the extraordinary case where the 
frequency ratio Vk :Ve :V, = 3 : 2 : 1 arises at the same radius -V3:2:i = 2.395. We can simply calculate that at 
;c= 2.369 the frequency ratio is VK:Ve:Vr = 3.05:2.019: 1 and atx=2.422, there is VK:V9:Vr = 2.95 : 1.981 : 1. 



a = 0.965729 a = 0.996383 




Fig. 3. Calculation of the width of the special value of the black hole spin parameter a3:2:i for which 
the multiple values frequencies should be observed. Black lines represent the ratio n.m = VK:Vr, VeiVr, 
dashed lines illustrate the relative frequency detuning S = 5% (i.e., n:m = 2±0.05 and 3 ±0.05). Left 
panel: The minimal value of the black hole spin parameter allowing the strong resonant phenomena where 
Vk : Ve :v, = 2.95 : 2.05 : 1 at a: = 2.556. Right panel: The maximal value of the black hole spin allowing the 
strong resonant phenomena where VK:Ve:Vj = 3.05 : 1.95 : 1 at a'= 2.255. 



and again the related four frequency ratio sets 

VK:ve: (VK-Vr) : (Vr = Ve-Vi) = 5:4:3:2, 

Vk : Ve :(Vr = Ve- v..) :(VK-Ve) =5:4:2:1, (52) 

and one five frequency ratio set 

VK:ve: (vk-Vi): (Vi =ve-Vr): (VK-Ve) = 5:4:3:2: 1. (53) 

Considering the extraordinary resonant spin a = fl3:2:i = 0.983, we can conclude 
that the resonant phenomena are possible in some region around the resonant radius 
•"^3:2:1 = 2.395 due to the allowed frequency detuning with detuning parameter o given 
by the frequency-response equations (see, e.g., Eq. |9] Eq. [T2b. The magnitude of the 
frequency scatter A = ea is fully determined by a concrete resonance under conside- 
ration and by physical conditions in the oscillatory system, nevertheless, we can put a 
general restriction on the allowed frequency detuning as few percent of the oscillator 
eigenfrequency (Nayfeh and Mook 1979) and use it to find region of the disc around 
•^3:2:1, where the resonance could appear We introduce the relations 

Vi=2vo±A, V2 = 3vo±A, (54) 

t Vl 

- = — =2±5, (55) 

u Vn 
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- = — = 3±5, (56) 

U Vo 

where 

5^-; (57) 
Vo 

the resonance could appear where the regions given by conditions ( fSSl l. ( |56]) overlap 
(see Fig.|2]i. We find that the maximal frequency detuning of 0.05 implies possibility of 
the strong resonant phenomena in the interval of xG (2.369;2.422) extended around the 
radius xt,-2:i — 2.395. Using this approach, we can give, for a given relative frequency 
detuning 5 = A/vq, the range of the black hole spin allowing the strong resonant phe- 
nomena. We simply find the spin where the overlap of the radial regions allowing 
resonance disappears (see Fig. |3]l. We choose characteristic values of 5 = 0.01,0.02 
and 0.05 and using a simple numerical code we determine the scatter around the spe- 
cial resonant spin values fl3:2:i, fl4:3:i and fl5:3:i. The results are summarized in Table[T] 
For the other special resonant spins and any appropriately chosen value of 5, the spin 
scatter can be given by the numerical code. (We give the spin with precision of 0.001 
as the 1% frequency scatter implies the spin scatter « 0.001.) 



Table 1 

Estimates of the spin value range where the strong resonant phenomena are possible for 
characteristic values of relative frequency detuning 8. 



8 [%] 


s:t:u = 3:2:l 


i:r:i( = 4:3: 1 


i- : f : w = 5 : 3 : 1 


1 
2 
5 


03:2:1 =U.yB3_oo03 
03:2:1 = 0.983_o (,05 
«3:2:1 = U yo3_Q q,, 


"4:3:1 = U-»DO_(, qq^ 

n sf;^;+o oio 

£14:3:1 =''■000-0.010 

n sfi^;+0 023 

04:3:1 ="'^°"-0.021 


n Qfio+O.OOl 

05:3:1 =0.962_oooi 

n Q(;o+0.003 

05:3:1 =0.962_ooq3 

n Q(;o+0 006 

05:3:1 =^-9b2_QQQj 



4 A Possible Application to Sgr A* 

The resonant phenomena could occur frequently, when different versions of the re- 
sonance model can be realized at a shared radius (or its close vicinity) fixed by the 
frequency ratio Vk : Vg : Vr = s : f : m. 

As explored above, there is an extraordinary resonant black hole spin fl3:2:i =0.983, 
when at the radius jc3:2:1 — 2.395 the frequency ratio is Vk : Ve : Vr = 3 : 2 : 1. Clearly, 
in vicinity of black holes with a = 0.983, the resonant phenomena should be strong, 
as the order of the resonances is of the lowest possible values, and, moreover, all the 
resonances, including those with beat frequencies, could cooperate efficiently even for 
frequencies slightly scattered from the exact resonant eigenfrequencies. The scatter of 
resonant frequencies strongly depends on the order of resonances and can be wide for 
resonances of very low order (Landau and Lifshitz 1976). 

When the simple combinational frequencies could be considered, the frequency 
ratios, e.g.,Vr:' (Vk — Ve) = 1:1, then allow strongest possible resonances. 

That is the reason why we could expect well observable QPOs in the field of black 
holes with spin close to the values allowing the cooperating resonant phenomena shar- 
ing a fixed radius. It is important to look for some candidate systems exploring the 
simple triple frequency ratio sets. In such situations we can determine the black hole 
spin with high precision, given by the frequency measurement error (see Fig.|4|. This 
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Fig. 4. En'or in determining the extraordinary resonant spin a3:2:i =0.983: tlie scatter of the black hole spin 
related to the 2% error in frequency measurements. It has to be confronted with the spin scatter allowing 
occuiTence of resonances (see Eqs (55) - (57)). Left panel: The case of top identity Vi,pp = Vk: the interval of 
allowed values of the spin is then a e (0.973; 0.992). Right panel: The case of bottom identity Vjown = Vr = 
(Ve - Vr): the interval of allowed values of the spin is a e (0.981; 0.984). 



could help very much in determining the other physical characteristics of the systems. 
Notice that errors of frequency measurements imply some errors in the spin determi- 
nation, as illustrated in Fig.|4] It depends also on the concrete resonances occurring at 
a given radius. Clearly, if the relevant frequency curves cross in a large (small) relative 
angle, the spin is determined with high (low) precision. 

The Galaxy centre source SgrA* can serve as a proper candidate system, since 
three QPOs were reported (but not fully accepted by the astrophysical community) for 
the system (Aschenbach 2004, Torok 2005) with frequency ratio corresponding to the 
extraordinary resonant spin 

(1/692) : (1/1130): (1/2178) w3: 2: 1 (58) 

and with the upper frequency being observed with a rather high error 

Vupp== (1.445 ±0.16) mHz. (59) 

Considering the standard epicyclic resonance model for the two upper frequencies 
Ve : Vr = 3 : 2, and the bottom frequency Vdown = Ve — Vr, we obtain Ve : Vr : (ve — Vr) = 
3:2: 1. Then identifying Vupp = Ve, we can show that for a < 1, the black hole mass 
M<2.3 X 10^ Mq, which is in clear disagreement with the allowed range of the SgrA* 
mass coming from the most recent analysis of the orbits of stars moving within 0.01 pc 
of SgrA* (Ghez et al. 2008) 

3.5xlO^M0<M<4.7xlO^M0 (60) 

and with the estimate M w 4 x 10^ obtained by Reid (2008). 

If the black hole is assumed to be at rest with respect to the Galaxy (/.e., has no mas- 
sive companion to induce motion), the fit could be further constrained to the interval 
(Ghez et al. 2008) 

4.1 X lO''M0<M<4.9x IO^'Mq (61) 

that is the most likely mass estimate at present, considering the discovery that SgrA* 
is nearly stationary at the Galactic center (Reid 2008). 

However, assuming a black hole with the spin comparable to the extraordinary 
resonant value of a sa 0.983, with the frequency ratio Vk : Ve : Vr = 3 : 2 : 1 at the sharing 
radius jC3:2:i =2.395, and identifying Vupp=VK, we obtain the black hole mass of SgrA* 
in the interval 

4.29X lO*'M0<M<5.36x IO'^Mq (62) 

which meets the allowed black hole mass interval (Eq.l60b at its high mass end and is 
in good agreement with the estimate (Eq.lMT). 
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5 Conclusions 

We present the possibilities for strong resonant phenomena arising in a shared radius 
in the field of black holes with spin appropriately tuned. We assume that the strong 
resonances could occur and influence each other causally, if Vk : Ve : Vr is in ratio of 
small integers s:t:u with s<5, when the order of the resonances could be low enough 
in order to enable strong resonance with relatively wide resonance frequency width 
(Landau and Lifshitz 1976). The strong resonances with lowest values of the frequency 
ratio occur for the extraordinary resonant spin a3:2:i = 0.983. 

Due to the allowed frequency detuning in the resonant phenomena (with maximal 
relative frequency detuning limited to few percent), the strong resonant phenomena 
could appear in black holes with spin values concentrated around the special resonant 
spin values. For a given relative frequency detuning 5, the spin scatter differs for dif- 
ferent special resonant spin values, being dependent on the position of the resonant 
radius (and the gradient of the frequency profiles at this radius), see Table [T] In the 
most interesting case of the extraordinary resonant spin fl3:2:i = 0.983, the spin scatter 
« ±0.003 for detuning 5 — 0.01 gives probably a realistic interval of black hole spins 
where the strong resonant phenomena could be observable. Of course, the value of 5 
(and the spin value scatter) depend on physical conditions and resonant phenomena in 
concrete sources and have to be discussed carefully for each individual source. 

There is an indication that the QPOs data observed in the Galactic centre source 
Sgr A* imply the black hole spin close to the extraordinary resonant value of aw 0.983, 
when Vk : Ve : Vr = 3 : 2 : 1 holds at the radius X3-2:i =2.395 and the black hole mass can 
be estimated to be in the interval of (4.8 ±0.5) x 10* M©. Therefore, it is interesting 
to check concordance of the QPOs observation induced data with a variety of the other 
observations of the Sgr A*. It should be stressed that more precise measurement of 
the QPOs frequencies will enable more precise determination of the black hole mass. 
Moreover, we expect NGC 5408 X-1 to be another candidate for a black hole admitting 
strong resonant phenomena because of the observed frequency set (Strohmayer et al. 
2007). 

It is important from the principal reason to note that if a black hole with appropri- 
ately tuned spin a admits strong resonant phenomena with different resonances sharing 
the same radius of the accretion disc, we could expect observable QPOs at resonant 
frequencies even in situations when the internal disc conditions are not convenient 
for starting up resonant phenomena, as observed, e.g., in Sgr A* source (Genzel et al. 
2003). 

On the other hand, the four (or five) QPOs frequencies observed in the microquasar 
GRS 1915H-105 can not be explained by the strong resonances model because of spe- 
cific distribution of the observed frequency ratios (two different pairs of frequencies 
with 3 : 2 ratio). In this case, the extended resonant model including the so called 
humpy oscillations induced by the "humpy" velocity profile of the accretion disc, as 
measured in the privileged family of locally non-rotating frames, can explain all the 
observed frequencies, if the black hole parameters are fixed at the spin a w 0.9998 and 
the mass M w 14.4 Mq (Stuchlik et al. 2007c), in agreement with the spectral fit of the 
spin a « 1 (McChntock et al. 2006). 

The conditions for strong resonant phenomena could be realized only for black 
holes with high values of dimensionless spin (a > 0.75). Therefore, the idea of strong 
resonant phenomena probably could not be extended to the neutron star systems, where 
we expect spin a <0.5, at least if the Hartle-Thorne metric parameters corresponding to 
the spin a and the quadrupole momentum q are close to the quasi-Kerr values with q ~ 
j^. Note that for the neutron star systems probably another version of the multiresonant 
idea could be realized in close agreement with observational data, where the so called 
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total precession resonance model with resonant frequencies Vk and V9 — Vr (related to 
the relativistic precession model, Stella and Vietri 1998) is involved and realized at 
different radii, as shown in Stuchlik et al. (2007e). 

In the field of Kerr black holes with a wO.866, the four frequency set ratio 4:3:2:1 
could be observed, if the resonant conditions allow combinational frequencies to be 
observable. For the spin flwO.962 or a wO. 882, even the five frequency set 5 : 4 : 3 : 2 : 1 
could be observable in some special circumstances. Clearly, it is not necessary that 
all the resonances are realized simultaneously and that the full five frequency set is 
observed at the same time. 

Undoubtly, it is worth to search for such special cases of frequency sets in observa- 
tional data as these could lead to precise determination of the black hole spin enabling a 
deeper understanding of a variety of astrophysical phenomena involved in determining 
behaviour of accretion discs in strong gravitational fields. We should note that when 
simple combinational frequencies not exceeding the Keplerian frequency are allowed, 
the lowest triple frequency ratio set 3 : 2 : 1 can be realized for black holes with the ex- 
traordinary resonant spin a w 0.983, but also for the black holes with three other values 
of the spin a = 0.866,0.882,0.962, if the uppermost frequencies are not observed for 
some reasons. 

Finally, it has to be stressed that observation of some characteristic frequency sets 
with low integer ratio could be a signature of a specific value of the black hole spin a 
enabling strong resonances at a fixed radius. However, generally, there exist few values 
of the spin a and the corresponding shared resonance radius allowed for a given fre- 
quency ratio set (Stuchlik et al. 2007b). Therefore, detailed analysis of the resonance 
phenomena, including the resonance strength and resonant frequency width, has to be 
considered in a concrete candidate system, and it has to be further confronted with the 
spin estimates coming from spectral analysis of the black hole system as given, e.g.,m 
McCHntock et al. (2006) and Middleton et al. (2006) for GRS 1915H-105, and Shafee 
et al. (2006) for GRO J1655-40, the line profiles (Fabian and Miniutti 2005, Dovciak et 
al. 2004, Zakharov 2003, Zakharov and Repin 2006), and the orbital periastron preces- 
sion of some stars moving in the region of Sgr A* (Kraniotis 2005, 2007), in order to 
establish the black hole spin. Very promising from this point of view seem to be studies 
of the energy dependencies of high-frequency QPO determining the QPO spectra at the 
QPO radii (Zycki et al. 2007) as they could be credited by the strong resonance condi- 
tions giving precisely the black hole spin and the radius where the observed QPOs are 
generated. 
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